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Abstract -Knowing the strategy of an opponent in a competitive environment conveys obvious 
evolutionary advantages. But this information is costly, and the benefit of being informed may not 
necessarily offset the additional cost. Here we introduce social dilemmas with informed strategies, 
and we show that this gives rise to two cyclically dominant triplets that form defensive alliances. 
The stability of these two alliances is determined by the rotation velocity of the strategies within 
each triplet. A weaker strategy in a faster rotating triplet can thus overcome an individually 
stronger competitor. Fascinating spatial patterns favor the dominance of a single defensive al¬ 
liance, but enable also the stable coexistence of both defensive alliances in very narrow regions of 
the parameter space. A continuous reentrant phase transition reveals before unseen complexity 
behind the stability of strategic alliances in evolutionary social dilemmas. 


Successful evolution of cooperation in the realm of so¬ 
cial dilemmas is a grand challenge that continues to at¬ 
tract research across the social and natural sciences [Ml- 
When individuals are torn between what is best for them 
and what is best for their society, cooperation can quickly 
come to play second fiddle to the pursuit of short-term 
personal benefits. Although studies in evolutionary game 
theory [ZHU] have revealed fundamental rules that pro¬ 
mote cooperation m, the chasm behind the Darwinian 
“only the fittest survive” and the abundance of coopera¬ 
tion in human and animal societies remains overwhelming. 

The availability of information has been the holy grail in 
modern game theory since its inception [13] . If one knew 
the rules of the game and what strategy the opponent will 
play, then the best response would be guaranteed [T4HT6I . 
In repeated evolutionary settings with many players and 
under the reasonable assumption that information about 
the opponents is not free, however, the prospects are much 
less clear. In fact, recent human experiments have revealed 
that, contrary to expectations, costly information about 
the neighbors can impair the evolution of cooperation in 


social dilemmas m- Important unanswered questions are 
thus: (i) Under which conditions are informed strategies 
players that invest into learning the strategies of other 
players — evolutionary stable? (ii) What are the proper¬ 
ties of these stable solutions? and (iii) Which mechanisms 
are responsible for their stability? 

In this letter, we answer these questions by introducing 
social dilemmas with informed strategies. The following 
matrix describes the payoffs of the four competing 
strategies: 
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Between unconditional cooperators (C) and defectors (D), 
in agreement with the traditional formulation of a social 
dilemma [18] . mutual cooperation yields the reward R, 
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Fig. 1: The direction of invasion between the competing strate¬ 
gies (see arrows), where the two cyclically dominant triplets 
Ic —► D —y C —y Ic and Ic —> Id —» C —y Ic form a defensive 
alliance against the external strategy (either Id or D). 


mutual defection leads to punishment P, and the mixed 
choice gives the cooperator the sucker’s payoff S and the 
defector the temptation T. Informed cooperators (Ic) and 
defectors (Id), on the other hand, invest e into knowing 
the strategy of their opponent and act accordingly. In 
particular, Ic players refuse to be exploited by defectors, 
while Id players avoid being punished when encountering 
other defectors. We suppose that being informed conveys 
an advantage to Id players compared to D players. With¬ 
out increasing the number of parameters further, the sim¬ 
plest way to ensure this is that D s are still being punished 
if they meet with an Id player. For simplicity but with¬ 
out loss of generality, we may decrease the number of free 
parameters by introducing T = 1+7", R = 1, P = —p 
and S = — r. Furthermore, the condition 0 < p < r, 
namely that the punishment for mutual defection is less 
than the sucker’s payoff and that the reward for mutual co¬ 
operation is smaller than the temptation to defect, yields 
T > R > P > S and thus returns the most demanding and 
widely studied social dilemma — the prisoner’s dilemma 
game (cf. [T9l - f32] ). 

With this definition of the payoffs, the directions of in¬ 
vasion between the competing strategies are as shown in 
Fig- [Q An unconditional cooperator (C) will always beat 
an informed cooperator (Ic), simply because the latter 
has to bear the cost e. True to the formulation of a tradi¬ 
tional social dilemma, an unconditional defector (D) will 
also beat C. This relation is reversed against Ic players 
due to their ability to avoid exploitation, as long as e < p. 
The Ic strategy is also superior to informed defectors (Id) 
because although both have to bear e, Ic players do bene¬ 
fit from mutual cooperation. Moreover, if we assume e < p 
and thus also e < r (because r > p), Id beats both D and 
C, respectively. 

The relations between the strategies shown in Fig. [Ore- 
veal the competition between two strategy triplets, namely 
Ic —^ D —y C —^ Ic and Ic —^ Id —^ C ^ Ic- These two 
triplets form a so-called defensive alliance [33H35] that pro¬ 
tects them against the external fourth strategy that is not 
within the loop. For example, the (Ic + Id + C) triplet 
is able to prevent strategy D from invading because both 


Ic and Id beat D. Similarly, the (Ic + D + C) prevents 
the invasion of Id- Despite of these relatively simple and 
clear-cut relations, several open questions remain, such as 
which informed strategy is able to prevail and how high 
can the cost e be, which defensive alliance is stable, or can 
the two triplets coexist? 

Since pattern formation always plays a decisive role in 
structured populations [36], we proceed with an in-depth 
analysis of the proposed evolutionary game on a L x L 
square lattice with periodic boundary conditions. Unless 
stated differently, we use random initial conditions such 
that all four strategies are uniformly distributed across 
the lattice. We simulate the evolutionary process in ac¬ 
cordance with the standard Monte Carlo simulation pro¬ 
cedure comprising the following elementary steps. First, 
a randomly selected player x acquires its payoff H x by 
playing the game with all its four neighbors. Next, player 
x randomly chooses one neighbor y 1 who then also ac¬ 
quires its payoff H y in the same way as previously player 
x. Lastly, player x adopts the strategy s y from player 
y with a probability determined by the Fermi function 
W(s y -y s x ) = (1 + exp[(n a; - II y )/AT]) -1 , where K = 0.1 
quantifies the uncertainty related to the strategy adop¬ 
tion process [Tj. The selected value ensures that better¬ 
performing players are readily followed by their neigh¬ 
bors, although adopting the strategy of a player that per¬ 
forms worse is not impossible either. Note that the results 
remain qualitatively similar in a wide range within the 
K < 2 interval. At very high K values, however, the 
system arrives to the random imitation limit. Each full 
Monte Carlo step (MCS) gives a chance for every player to 
change its strategy once on average. During simulations, 
we have used lattices with up to 4000 x 4000 players to 
avoid undesired finite-size effects, and we have determined 
the stationary fractions of the strategies in the stationary 
state after sufficiently long relaxation times extending up 
to 4 • 10 5 MCS. In addition, we have repeated every run 
up to 100 times with different initial conditions to further 
decrease statistical fluctuations. 

In Fig. [2] we show the complete r — p phase diagram 
of the system as obtained for a representative value of e. 
We would like to stress, however, that qualitatively similar 
behavior can be obtained for all relevant values of param¬ 
eters e or K. Several interesting solutions deserve to be 
highlighted. First, the Ic strategy is always viable, even 
though in its absence defectors would easily eradicate un¬ 
conditional cooperators. However, the later also become 
viable as second-order free-riders of Ic (first-order free¬ 
riders are either D or Id players). For low values of p 
the (Ic + D + C) triplet prevails and Id players die out, 
while for high values of p the (Ic + Id + C) triplet wins 
on the expense of D. Interestingly, there is a very nar¬ 
row dividing stripe where both triplets coexist, and where 
thus all four strategies survive. The dependence on r is 
even more interesting. At an intermediate value of p , the 
system undergoes a succession of four consecutive phase 
transitions as r increases, after which it ends in the same 
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Fig. 2: Full r — p phase diagram, as obtained for e = 0.3. 
Solid red line denotes continuous phase transitions from the 
very narrow but stable Ic + D + C + Id phase to the stable 
Ic + D + C phase, while the dotted blue line denotes reentrant 
continuous phase transitions to the stable Ic + Id + C phase. 
Dashed green line is the p < r border. 


Ic+Id+C phase in which it started. This reentrant phase 
transition is illustrated in the top panel of Fig. [3j where 
both Ic + Id + C phases are separated by two very narrow 
Ic + D + C + Id stripes and a broader Ic + D + C region 
in the middle. What is truly remarkable is that the com¬ 
petition between D and Id (which exchange in the two 
triplets) is not directly affected by the value of r. What 
is more, four out of all the six relations in the food web 
are independent of r. While a reentrant phase transition 
has been reported before in the realm of predator-prey 
interactions due to rewiring 132 !, here it is due solely to 
strategic interactions. 

An explanation for the unexpected behavior thus has to 
lie in the D —> C and Id —> C invasion rates, which are 
directly affected by r. The bottom panel of Fig. [3] reveals 
a monotonous dependence of the two invasion rates on r, 
whereby both curves start and terminate at an identical 
level. Thus, at both ends the two triplets have an equally 
fast internal dynamics, but since Id players beat D play¬ 
ers, the dominance belongs to the (Ic + Id + C) triplet. 
But for intermediate values of r, the D —> C invasion is 
significantly more efficient than the Id ~+ C invasion, and 
the average rotation of strategies within the (Ic + D + C) 
triplet is therefore much faster than in the (Ic + Id + C) 
triplet. Due to this difference in the internal dynamics, the 
Id players meet with their direct predator, the Ic players 
(see the food web in Fig. [TJ, more frequently. Because of 
this the faster rotating alliance is able to overcome the di- 




Fig. 3: Top panel shows a characteristic cross-section of the 
phase diagram in Fig. [3] as obtained for p = 0.45. Depicted 
are stationary fractions of the four competing strategies in de¬ 
pendence on r. After a succession of four continuous phase 
transitions, the system ends up in the same Ic + Id + C phase 
in which it started. The reentrant phase transition on both 
sides is preceded by a very narrow region (marked by arrows) 
where all four strategies coexist. Bottom panel shows the inva¬ 
sion rates between cooperator and defector strategies (see leg¬ 
end) in dependence on r, as obtained for e = 0.3 and p = 0.45. 
The hysteresis is crucial for the understanding of the reentrant 
phase transition to the Ic + Id+C phase (see main text for 
details). 


rect evolutionary disadvantage that D players have against 
Id players, and in fact reverse the outcome of the compe¬ 
tition between the two defensive alliances in favor of the 
(Ic + D + C) triplet. Naturally, as soon as the two inva¬ 
sion rates meet at either the low or the high limit of r, the 
difference in the rotation speeds within the alliances dis¬ 
appears, and the original victor, the (Ic + Id+C ) triplet, 
is restored. 

As we have shown, when the D —>• C and Id ~+ C 
invasion rates are equal, the winner is decided by virtue 
of the direct dominance of Id players over D players and 
thus the (Ic + Id + C) triplet wins. On the other hand, 
when the D —> C invasion is significantly more efficient, 
then the faster average rotation of strategies within the 
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Fig. 4: Consecutive snapshots of the square lattice, as obtained for r = 0.65 (top row) and r = 0.95 (bottom row) from a 
prepared initial state. During the relaxation period, each half of the lattice contained the strategies of one triplet only. After 
the characteristic patterns evolved (leftmost panels), the separating wall (vertical white line) was removed and the two defensive 
alliances could then start to compete with each other. In the top row, the (Ic + D + C) triplet, which is depicted by dark 
blue, light blue, and light red colors (overall lighter), gradually compresses the ( Ic + Id + C) triplet, which is depicted by dark 
blue, light blue, and dark red colors (overall darker). For convenience, the typical size of the shrinking alliance is encircled 
with a white ellipse. The snapshots were taken at t = 0, 100, 200 and 500 MCS from left to right. In the bottom row, the 
(Ic + Id + C) triplet gradually rises to dominance because Id players are superior to D players. Accordingly, the domains of 
the (Ic + D + C) triplet (encircled with white ellipses for clarity) disappear over time. Here the snapshots were taken at t = 0, 
100, 1000 and 4500 MCS. Note that the evolutionary process in this case is almost ten times slower than in the top row. Other 
parameter values used in both cases are L = 200, p = 0.4 and e = 0.3. 


(Ic + D + C) triplet proves fatal for the (Ic + Id + C) 
triplet. These two different mechanisms that decide the 
competition between the two defensive alliances can also 
be appreciated through patterns that form on the square 
lattice, as illustrated in Fig. [I] For clarity, we have split 
the square lattice in half and used prepared initial states 
to highlight the competition between the two triplets. 

The top row of Fig. [I] depicts the evolution when the 
(Ic + D + C) is the more dynamic triplet. It can be ob¬ 
served that the other triplet gradually decreases in size 
because its defensive capacity is compromised. To illus¬ 
trate this monotonous decay, we have marked the typical 
size of the shrinking area by white ellipses during the in¬ 
termediate stages of the evolutionary process. Although 
we only present a few snapshots, the difference in the ro¬ 
tation velocity within the two defensive alliances is still 
remarkably obvious. While the large majority of the typ¬ 
ical pattern that characterizes the (Ic + Id + C) triplet 
(right side of the lattice) remains practically unchanged, 
the (Ic + D + C) triplet (left side of the lattice) com¬ 
pletes a whole strategy rotation during this time. Note 


that in the leftmost panel the majority of players within 
the (Ic + D + C) alliance is in the cooperator state, in 
the second panel D players abound, and then in the third 
panel the original fraction of C + Ic players is again re¬ 
stored. Actually, such an “oscillation” is just a finite-size 
effect that is very useful for illustrative purposes, while in 
fact the fractions of the strategies are time-independent 
in the stationary state if a sufficiently large system size is 
used. 

Conversely, the bottom row of Fig. [4] depicts the evo¬ 
lution when the strategies within the two triplets rotate 
equally fast. Accordingly, the pattern formation is sig¬ 
nificantly different from the one shown in the top row 
of Fig. 01 Here the roles of strategies Id and D in the 
two alliances are interchangeable — both strategies beat 
strategy C with practically the same efficiency. As a con¬ 
sequence, the patterns formed by the two triplets are very 
similar, and this not just approximately but quantitatively 
as well. If we compare the fractions of D and Id players 
in the stable Ic + D + C and Ic + Id + C phases around 
the right reentrant phase transition depicted in the top 
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Fig. 5: Top panel shows the fraction of strategy D in the Ic + 
D + C phase at p = 0.33 and the fraction of the Id strategy in 
the Ic+Id+C phase at p = 0.6 in dependence on r, as obtained 
for e = 0.3. The bottom panel shows the corresponding fraction 
of Ic players {pi c ) in the two triplets at the same two values 
of p (see legend). Counterintuitively, due to cyclic dominance, 
as the temptation to defect increases, the fraction of defectors 
goes down while the fraction of informed cooperators goes up. 


panel of Fig. O we can observe that they are practically 
equivalent (just switch values). Indeed, the superiority 
of the Id strategy manifests only when it is faced with 
strategy D. Because of this “individual” competition, the 
evolutionary process is much slower than in the top row 
(note that it takes ten times longer to reach the single- 
triple stationary state). Another more subtle consequence 
of identically fast average rotations of strategies within the 
two competing triplets is that there is a good chance for 
a lucky D player to emerge and spread efficiently on the 
expense of C players. In doing so, regions of the square 
lattice that already belonged to the {Ic + Id + C) triplet 
may be temporarily lost to the {Ic + D + C) triplet. We 
have marked such unexpected patches of the {Ic + D + C ) 
alliance by white ellipses in the third panel of the bottom 
row of Fig. a 

An unexpected but important consequence of cyclic 
dominance in the studied system is also the peculiar de¬ 
pendence of the fractions of individual strategies within 


the two defensive alliances on r, which we illustrate in 
Fig. [5] Based on the payoff matrix, we would expect 
that increasing values of r will unambiguously support 
both defector strategies and at the same time weaken un¬ 
conditional cooperators, as is the case in the traditional 
two-strategy prisoner’s dilemma game. But as results pre¬ 
sented in Fig.[5]show, exactly the opposite happens regard¬ 
less of which defensive alliance is dominant. This coun¬ 
terintuitive behavior is a beautiful consequence of cyclic 
dominance, where the direct support for a particular strat¬ 
egy will actually promote its “superior” in the food web. 
Such evolutionary outcomes are common in rock-paper- 
scissors-like systems, where feeding a prey often promotes 
the predator [38H45] . In the studied social dilemma with 
informed strategies this is exactly what happens, since in¬ 
creasing the value of r will directly support D and Id 
strategies, which however are both prey to the Ic strat¬ 
egy (see the food web in Fig.[T|). As demonstrated in the 
bottom panel of Fig. 0 indeed the fraction of Ic goes up 
in both phases as r increases. An additional consequence 
of increasing the value of r is that it lowers the fitness of C 
players. Since the latter act as second-order free-riders of 
Ic players, and are at the same time prey to both D and 
Id strategies, this also promotes the viability of informed 
cooperators. 

Before summarizing our findings, we briefly discuss the 
role of the information cost e. Throughout this work, we 
have considered e = 0.3, as it yields the most representa¬ 
tive system dynamics. By increasing the value of e, the 
informed strategies become weaker. While the Id —> D re¬ 
lation still holds, the intensity of this invasion decreases. 
Consequently, the winning mechanism of the {Ic+Id+C) 
triplet is negatively affected, and the region where the 
faster rotating {Ic + D + C) triplet is able to compensate 
the Id — > D relation expands. For e > 0.64 the Ic+Id+C 
alliance disappears completely. On the other hand, lower e 
values intensify the Id —> D invasion, and accordingly, the 
region where the Ic + Id + C triplet is dominant expands 
on the expense of the Ic + D + C triplet. Nevertheless, the 
faster rotating Ic + D + C alliance can defend its stability 
throughout the small e region. 

In summary, we have proposed and studied an evolu¬ 
tionary social dilemma game with informed strategies. We 
have identified elementary relations between the four com¬ 
peting strategies in the governing food web, which revealed 
the existence of two three-strategy defensive alliances. By 
means of systematic Monte Carlo simulations, we have re¬ 
vealed key determinants of stability of the two competing 
triplets, and we have observed fascinating solutions that 
are driven by pattern formation and cyclic dominance. 
We have shown that a direct evolutionary advantage of 
a strategy within a defensive alliance can be compensated 
by the other alliance through a faster internal rotation of 
its strategies. Thus, even though in the food web the in¬ 
formed defectors are superior to unconditional defectors, 
the alliance whose defense relies on the weaker strategy 
can still prevail. In particular, we have shown that the 
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competition between a direct food-web-based evolutionary 
advantage and an evolutionary advantage that is rooted in 
the spatiotemporal dynamics of a defensive alliance gives 
rise to a reentrant phase transition. In-between the two 
stable three-strategy phases that form the reentrant pair, 
we have also identified a very narrow region of coexis¬ 
tence of both defensive alliances, which emerges as a con¬ 
sequence of a delicate equilibrium of the two competing 
mechanisms. 

When approaching more realistic conditions to study 
the evolution of cooperation [46] [47], the emergence of 
cyclic dominance becomes more and more likely. Our re¬ 
search highlights that the final outcome in such cases de¬ 
pends sensitively not just on the individual relations be¬ 
tween the competitors as determined by payoff elements, 
but in an equally strong manner also on the dynamical 
properties of alliances. This in turn makes interventions 
aimed as steering such complex systems risky and to a 
large degree unpredictable. 

* * * 

This research was supported by the Hungarian National 
Research Fund (Grant K-101490), the Slovenian Research 
Agency (Grant P5-0027), and the Deanship of Scientific 
Research, King Abdulaziz University (Grant 76-130-35- 
HiCi). 

REFERENCES 

[1] Szabo G. and Fath G., Phys. Rep., 446 (2007) 97. 

[2] Santos F. C., Pinheiro F., Lenaerts T. and Pacheco 
J. M., J. Theor. Biol. , 299 (2012) 88. 

[3] Nowak M. A., J. Theor. Biol., 299 (2012) 1. 

[4] Perc M. and Szolnoki A., BioSystems, 99 (2010) 109. 

[5] Rand D. A. and Nowak M. A., Trends in Cognitive 
Sciences, 17 (2013) 413. 

[6] Pacheco J. M., Vasconcelos V. V. and Santos F. C., 
Physics of Life Reviews, 11 (2014) 573. 

[7] Mesterton-Gibbons M., An Introduction to Game- 
Theoretic Modelling, 2nd Edition (American Mathemati¬ 
cal Society, Providence, RI) 2001. 

[8] Nowak M. A., Evolutionary Dynamics (Harvard Univer¬ 
sity Press, Cambridge, MA) 2006. 

[9] Sigmund K., The Calculus of Selfishness (Princeton Uni¬ 
versity Press, Princeton, NJ) 2010. 

[10] Lugo H. and San Miguel M., Sci. Rep., 5 (2015) 7776. 

[11] Wu Z.-X. and Yang H.-X., Phys. Rev. E, 89 (2014) 
012109. 

[12] Nowak M. A., Science, 314 (2006) 1560. 

[13] von Neumann J. and Morgenstern O., Theory of 
Games and Economic Behaviour (Princeton University 
Press, Princeton, NJ) 1944. 

[14] Macy M. W. and Flache A., Proc. Natl. Acad. Sci. 
USA, 99 (2002) 7229 

[15] Rustagi D., Engel S. and Kosfeld M., Science, 330 
(2010) 961 

[16] Szolnoki A. and Perc M., Phys. Rev. E, 85 (2012) 
026104 


[17] Antonioni A., Cacault M. P., Lalive R. and 
Tomassini M., PLoS ONE, 9 (2014) ell0788. 

[18] Santos F. C., Pacheco J. M. and Lenaerts T., Proc. 
Natl. Acad. Sci. USA, 103 (2006) 3490. 

[19] Zimmermann M. G., Egui'luz V. M. and San Miguel 
M., Phys. Rev. E, 69 (2004) 065102(R). 

[20] Santos F. C. and Pacheco J. M., Phys. Rev. Lett., 95 
(2005) 098104. 

[21] Gomez-Gardenes J., Campillo M., Flori'a L. M. and 
Moreno Y., Phys. Rev. Lett., 98 (2007) 108103. 

[22] Ohtsuki H., Nowak M. A. and Pacheco J. M., Phys. 
Rev. Lett., 98 (2007) 108106. 

[23] Fu F., Hauert C., Nowak M. A. and Wang L., Phys. 
Rev. E, 78 (2008) 026117. 

[24] Poncela J., Gomez-Gardenes J., Flori'a L. M., 
Moreno Y. and Sanchez A., EPL, 88 (2009) 38003. 

[25] Fu F., Wu T. and Wang L., Phys. Rev. E, 79 (2009) 
036101. 

[26] Fu F., Nowak M. A. and Hauert C., J. Theor. Biol, 
266 (2010) 358. 

[27] Jiang L.-L., Wang W.-X., Lai Y.-C. and Wang B.-H., 
Phys. Rev. E, 81 (2010) 036108. 

[28] Antonioni A. and Tomassini M., PLoS ONE, 6 (2011) 
e25555. 

[29] Lee S., Holme P. and Wu Z.-X., Phys. Rev. Lett., 106 
(2011) 028702. 

[30] Tanimoto J., Brede M. and Yamauchi A., Phys. Rev. 
E, 85 (2012) 032101. 

[31] Hilbe C., Nowak M. and Sigmund K., Proc. Natl. Acad. 
Sci. USA, 110 (2013) 6913. 

[32] Santos M., Dorogovtsev S. N. and Mendes J. F. F., 
Sci. Rep., 4 (2014) 4436. 

[33] Szabo G. and Czaran T., Phys. Rev. E, 64 (2001) 
042902. 

[34] Szabo G., J. Phys. A: Math. Gen., 38 (2005) 6689. 

[35] Kim B. J., Liu J., Um J. and Lee S.-L, Phys. Rev. E, 
72 (2005) 041906. 

[36] Szolnoki A., Mobilia M., Jiang L.-L., Szczesny B., 
Rucklidge A. M. and Perc M., J. R. Soc. Interface, 11 
(2014) 20140735. 

[37] Han S.-G., Park S.-C. and Kim B.-J., Phys. Rev. E, 79 
(2009) 066114. 

[38] Frean M. and Abraham E. D., Proc. R. Soc. Lond. B, 
268 (2001) 1323. 

[39] Dobrinevski A., Alava M., Reichenbacli T. and 
Frey E., Phys. Rev. E, 89 (2014) 012721. 

[40] Ni X., Wang W.-X., Lai Y.-C. and Grebogi C., Phys. 
Rev. E, 82 (2010) 066211. 

[41] Wang W.-X., Lai Y.-C. and Grebogi C., Phys. Rev. 
E, 81 (2010) 046113. 

[42] Rulands S., Zielinski A. and Frey E., Phys. Rev. E, 

87 (2013) 052710. 

[43] Vukov J., Szolnoki A. and Szabo G., Phys. Rev. E, 

88 (2013) 022123. 

[44] Laird R. A. and Schamp B. S., J. Theor. Biol, 365 
(2015) 149. 

[45] Groselj D., Jenko F. and Frey E., Phys. Rev. E, 91 
(2015) 033009. 

[46] Capraro V., Smyth C., Mylona K. and Niblo G. A., 
PLoS ONE, 9 (2014) el02881. 

[47] Hauser O. P., Rand D. G., Peysakhovich A. and 
Nowak M. A., Nature, 511 (2014) 220. 


P-6 



